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Survival probability of surface excitations in a 2d lattice: non-Markovian effects and 

Survival Collapse. 

E. Rufeil FiorO and H. M. PastawskQ 

Facultad de Matemdtica, Astronomia y Fisica, Universidad Nacional 

de Cordoba, Ciudad Universitaria, 5000 Cordoba, Argentina. 

The evolution of a surface excitation in a two dimentional model is analyzed. I) It starts quadrat- 
ically up to a spreading time ts- 11) It follows an exponential behavior governed by a self- consistent 
Fermi Colden Rule. Ill) At longer times, the exponential is overrun by an inverse power law de- 
scribing return processes governed by quantum diffusion. At this last transition time tn a survival 
collapse becomes possible, bringing the survival probability down by several orders of magnitude. 
We identify this strongly destructive interference as an antiresonance in the time domain. 



I. INTRODUCTION 



We consider the dynamics of a charge excitation in a 
typical model for Tamm states 1]. Similar tight-binding 
models |3 are used to describe a variety of situations: 
molecules absorbed in metallic substrates, decay of high- 
energy electron excitations and the decoherence caused 
by a weak interaction with an "environment" whose spec- 
trum is dense. The decay of the survival probability 
^00 (^) of the resulting resonant state is usually described, 
within a Markovian approximation, by the Fermi Golden 
Rule (FGR). However, this description contains approx- 
imations that leave aside some intrinsically quantum be- 
haviors. Various works on models for nuclei, composite 
particles 0: S H ^ excited atoms either in a free electro- 
magnetic field p or in photonic lattices 7] , showed that 
the exponential decay has superimposed beats and does 
not hold for very short and very long times, compared 
with the lifetime of the system. 

In Ref. 8] we presented a model describing the evo- 
lution of a surface excitation in a semi-infinite chain, a 
model that is solved analytically and susceptible for an 
experimental test J2| . Here, we present a general analysis 
showing the quantum nature of the deviations from the 
Fermi Golden Rule. Then, we numerically solve a model 
consisting of an excited add atom in a two dimensional 
lattice. We identify three time regimes in the decay of 
the survival probability Poo{t): 1) For short times the de- 
cay is quadratic, as is expected when the coupling of the 
local state with the continuum is perturbative. 2) An 
intermediate regime characterized by an exponential be- 
havior, the self- consistent Fermi Golden Rule (SC-FGR) 
where the rate, the pre-exponential factor and the char- 
acteristic frecuency are found self-consistently. 3) A long- 
time regime in which the exponential decay of the pure 
survival probability is overrun by an inverse power law, 
which is identified with the return probability enabled by 



the slow quantum diffusion in the substrate. At this last 
cross-over, the oscillations could lead to a dip in Poo{t) 
of several orders of magnitude. This survival collapse 
is identified with a destructive interference between the 
pure survival amplitude, i.e., the SC-FGR component, 
and the return amplitude, associated with high orders in 
a perturbation theory. 



II. SURVIVAL PROBABILITY OF A SURFACE 
EXCITATION 

We consider the evolution of a surface excitation, pre- 
pared in the state |0) in a Hamiltonian with finite spec- 
trum as is the case of most excitations in a lattice. The 
survival probability is 



Poo{t)= {0\exp[-iHt/h]\0)0{t) 
^h'\G^,{t)\\ 



where 



gS, (t) 



I^G^oie)exp[-iet/h], 



(1) 
(2) 

(3) 



is the retarded Green's function for a single fermion. 
Expanding the initial condition in the eigenstates of H 
yields |l|l3 



-Poo(t) 



e (t) / d£ A^o (e) exp[-iet/h] 



(4) 
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where Nq {e) is the Local Density of States (LDoS) at 
site 0th which, in terms of the retarded Green's function 
is Nq (e) = — 1 /Trim G^(6:). Hence, we can evaluate the 
survival probability using the Fourier transform of the 
LDoS, which can be accurately calculated in the energy 
representation where the integral is limited to the spec- 
tral support. Besides, a clear identification of quantum 
interferences will be obtained by analyzing the argument 
of the square modulus. 

In order to evaluate the local dynamics, we perform the 
integral in Eq. Q using the residue theorem and following 
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FIG. 1: Local spectrum (LDoS) in the complex plane z = e+i 
e. sl and £u are the lower and upper band-edges, respectively. 
The pole appears in £r — iFo. The integration path is shown 
with dotted lines; consist of four straight lines and two arcs, 
that avoid the band-edges singularities. 



the path shown in the Fig. ^ I^i the analytical continu- 
ation No{z) = A^o(^ + i^O^ resonances appear as poles in 
the complex plane. We will consider Hamiltonians where 
an initially unperturbed state of energy sq = /o|i^|0 ) 
interacting with a continuum is a well defined resonance, 
i.e., the expansion of |0) in terms of the eigenstates has 
a small breath Fq around an energy Sr = £0 + Aq, where 
Ao = A (e = Sr) is a small shift due to the interaction. 
We exclude systems with localized eigenstates [ll| . Then 
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(5) 



return correction from quantum diffusion 



where 



a = 27ri lim [{z - Sr ^ iTq) No{z)] , (6) 
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The first term of Eq.Q already supersedes the usual 
Fermi Golden Rule approximation since it has a pre- 
exponential factor (|ap > 1) and the exact rate of decay 
Fq. This result is the self- consistent Fermi Golden Rule 
(SC-FGR). By analogy with a classical Markov model, 
this exponential term is identified with a ''pure survivaF 
amplitude. Within the same analogy^ the second term 
will be called ''return^^ amplitude, as it is fed upon the ini- 
tial decay. The first term is the dominant one for a wide 



range of times, while the "quantum diffusion" described 
by the second, dominates for long times and brings out 
the details of the spectral structure of the system. The 
second term is also fundamental for the normalization at 
very short times where the most excited energy states of 
the whole system can be virtually explored. Both terms 
combine to provide the initial quadratic decay required 
by the perturbation theory: 
+2 



Poo (t) 



t^ 



((^-^^)% + 



(8) 



Here ((e — e^)^)^ is the energy second moment of the 
density Nq {e). This expansion holds for a time shorter 
than the spreading time ts of the wave packet formed 
upon decay. 

For long times, the behavior of Poo {t) is governed by 
the slowly decaying second term in Eq.©. Only small 
values of e^ contribute to the integral. This restricts 
the integration of the LDoS to a range near the band- 
edges. Then, one can go back to Eq.(^ and perform the 
Fourier transform retaining only the van Hove singulari- 
ties [1^ [13 at these edges. The relative participation of 
the energy states at each edge of the LDoS is given by 
the relative weight of the Lorentzian tails at these edges 
p= [{er-eL)'^^Tl]/[{eu-er)'^^rl]. Then, the survival 
probability for long times is 

Poo{t) ^ [1 + /?' - 2(3cos{Bt/h)] 



fde'e-'''^^No{eL-ie') 



(9) 



This means that the long time behavior is just the power 
law decay of the integral multiplied by a factor containing 
a modulation with frequency B/h. 



III. SURVIVAL COLLAPSE 

In steady state transport ^12] as well as in dynamical 
electron transfer [13 there are situations when a par- 
ticle can reach the final state following two alternative 
pathways. Since each of them collects a different phase, 
this allows a destructive interference blocking the final 
state. This phenomenon has been dubbed antiresonance 
[1^ \l^ . It extends the Fano resonances describing the 
anomalous ionization cross-section [l^. In the present 
case, the survival of the local excitation also recognizes 
two alternative pathways: the pure survival amplitude, 
which is typically described by the Fermi Golden Rule, 
and the pathways where the excitation has decayed, ex- 
plored the substrate, and then returns. These two al- 
ternatives can interfere. We rewrite Eq.© to emphasize 
that the survival probability Poo (t) is the result of two 
different contributions: 



Poo {t) 



= 1^ 



^5 + ^i? 

2 



*flr + 2Re[*J*fl], 



(10) 
(11) 




FIG. 2: 2d-lattice with an add atom with site energy eo and 
hopping Vb. 



where the phase in ^n arise from the exponentials with 
sl and eu (the LDoS is real for any argument). Hence, 



^-rot/h^-iis^-SL)t/h 



*s(i) = |a|e 

*ii(t) = |*flW|e'^(*); 

0ft) = arctan( ^fg"jf^^^^)^^ ). 



(12) 
(13) 

(14) 



where Eq. (|T^ results using the long time limit of Eq.Q. 
While the interference term in Poo {t) is present along the 
whole exponential regime, it becomes important when 
both, the pure survival amplitude and the return contri- 
bution, are of the same order. This occurs at the cross- 
over time tn between the exponential regime and the 
power law. The interference term can produce a survival 
collapse^ i.e., a pronounced dip that takes Pqo {t) close to 
zero (see Fig. |2|. In order to obtain a full collapse, two 
simultaneous conditions are needed; 

\^s{tR)\ = \^R{tR)\ and (15) 

{er - Sl) tn/h - (j) (tn) = (tt - (/)«) + 27rn, (16) 

which are satisfied with a fair precision because the re- 
turn amplitude has a phase with a slow variation: 

\{er - sl) /h\ > To/h > 2n/tR > \^ {tR)\ liR, (17) 

while, the pure survival term oscillates rapidly. When 
both amplitudes are of the same order, the destructive 
interference will be noticeable. 



IV. DECAY IN A 2D-SYSTEM 

The above results (Eq.(|S|), Eq.Q, and the survival 
collapse effect) were verified and quantified in a recent 
publication 8], which solves the dynamics of a surface 
spin excitation weakly coupled to a semi-infinite chain of 
interacting spins. The integral in Eq.Q is analytically 
solved for the different time regimes (short, exponential 
and long time). Also, the cross-over from the short time 
regime to the exponential SC-FGR t^, and the cross- 
over from the SC-FGR to the power law regime Ir^ were 
found. The survival collapse takes place at Ir. 



Here, we shall consider a square lattice with an add 
atom. The Tight Binding Hamiltonian is 



^ = X^ 1^) ^n (^1 - ^ 1^) K,m (^1 



(18) 



where each (r|n) is centered around the corresponding 
lattice site n, Sn are the site energies and Vn,m are the 
hoppings. We consider the case where the 0th site de- 
notes the add atom, i.e., is different from the others 
sites in both site energy, ^q 7^ ^n = 41/ and hopping 
Vb,i = Vb < Vn.m = y ' This defines a continuous spec- 
trum in the range [0, P = ^V\, The Green function of 
this problem, and hence the LDoS, is evaluated using 
the Dyson equation 



[Gol, (.)]"' 



gT (e) 



Vo,iGfl°He)Vi,o, 



(19) 



following the general continued fraction procedure de- 
scribed in Ref. 2]: 



G^oie) 



1 



^0 



Vo'GTiey 



(20) 



^i?(0). 



where Gii"" (s) is the Green function for a periodic square 
lattice (Ref. U^). 

For this system, the local second moment of the Hamil- 
tonian is Vq. The short time regime (Eq.Q) holds up to 
a time ts in which the quadratic decay becomes an ex- 
ponential. A good estimate of ts is obtained from the 
minimal distance between the short time decay and the 
exponential decay. We can use the first pole approxi- 
mation (evaluating Gii {e = eo) in Ea. (|20j) ) to obtain 
Fo ^ 7tV^n[^^ {eo) , which coincides with the FGR. In 
the same order of approximation we can take |ap ^ 1 in 
Eq.((S|). Then we obtain 



ts « hnN[°^ (£o) . 



(21) 



This result shows that the spreading time ts is only de- 
termined by N\^ ^ (^o) 5 the local density of states at the 
first site of the unperturbed substrate. 

We verify Eq.(|5|), Eq.Q, and the survival collapse 
effect, using the analytical expression for Gqo (^) ^^^ 
No{e), and performing the numerical Fourier transform. 
In Fig. |31we show Poo(^)- The curve shows the exponen- 
tial SC-FGR, which then is overrun by a 1/t^ power law 
decay. The cross-over time tR is easily identified through 
the survival collapse shown as a dip in Poo- There, the 
survival probability suddenly decreases from its average 
by almost three orders of magnitude. The inset shows 
the small oscillation that modulates the power law. 

It is important to note that the square power law de- 
cay obtained for long times is a consequence of the 6 {e) 
dependence of the LDoS (inset of Fig. Q), i.e., this power 
law is consistent with Eq.Q taken together with Eq.(|2n|). 
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FIG. 3: Local polarization, in a doble logarithmic scale, as 
a function of time. We consider an unperturbed energy of 
£o/V = 2 and interaction strength Vo/V = 0.4. This is the 
case that we consider in FigQ The decay exhibits: The ex- 
ponential behavior as described by the self- consistent Fermi 
Golden Rule, and an asymptotic square power law decay. The 
inset shows the oscillation that modulates this decay. The 
cross-over time tn when the survival collapse takes place is 
indicated. 



V. CONCLUSIONS 

In the present work we have discussed the dynamics 
of a local excitation that decays through a weak interac- 
tion with a continuum spectrum with finite support. Our 
approach goes beyond the usual Markovian approxima- 
tion that uses the Fermi Golden Rule to describe these 
environmental interactions. 

The evolution starts with the expected quadratic de- 
cay. Then, it follows the usual exponential FGR regime, 
but with a corrected rate and a pre-exponential factor, 
i.e., the SC-FGR. Finally, we get the long time regime, 
that consists of a square power law decay modulated by 
oscillations whose frequency is determined by the band- 
width. This power law decay is a consequence of the (e) 
behavior of the LDoS in the band-edge (Eq.Q) and is 
identified with the quantum diffusion in the substrate. 
Hence, anomalies in the excitation decay gives informa- 
tion about the substrate dynamics. 

Finally, we predict the existence of the survival col- 
lapse. This non-Markovian result fully considers the 
memory effects to infinite order. Such effect, hinted but 
not explained in previous works, is visualized as the de- 
structive interference between the pure survival ampli- 
tude and the return amplitude that arises from pathways 
that have already explored the rest of the system. 
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